A short-range correlation energy density functional with multi-determinantal 

reference 



Julicn Toulouse, Paola Gori-Giorgi, and Andreas Savir0 
Laboratoire de Chimie Theorique, 
CNRS et Universite Pierre et Mane Curie, 
4 place Jussieu, 75252 Paris, France. 
(Dated: February 2, 2008) 

We introduce a short-range correlation density functional defined with respect to a multi- 
determinantal reference which is meant to be used in a multi-determinantal extension of the Kohn- 
Sham scheme of density functional theory based on a long-range/short-range decomposition of the 
Coulomb electron-electron interaction. We construct the local density approximation for this func- 
tional and discuss its performance on the He atom. 



I. INTRODUCTION 

One of the main difficulties in the Kohn-Sham (KS) [1] 
scheme of density functional theory (DFT) 2] is to find 
approximations for the exchange-correlation energy func- 
tional that correctly describe (near-)degeneracy or long- 
range (e.g., van der Waals) correlation effects. To cir- 
cumvent this difficulty, a multi-determinantal extension 
of the KS scheme based on a long-range/short-range de- 
composition of the Coulomb electron-electron interaction 
has been proposed [3, i, i, 0, 0, i, i, [13, [H d • The 
idea behind this separation is that correlation effects due 
to the short-range part, involving the correlation cusp, 
could well be described by the local density approxima- 
tion (appropriately modified) ; correlation connected with 
the long-range part could well be dealt with using stan- 
dard wave-function methods of quantum chemistry. 

In this approach, the ground-state energy of a iV- 
electron system in a nuclei-electron potential f„e(r) is 
obtained in principle exactly by minimization over multi- 
determinantal wave functions ^' 

E = min((*|f + 14e + l^ie'^|*) 

+i?^^^''[n*] + + (1) 

where T is the kinetic energy operator, Vne — 
'^iVne{TCi) IS the uuclei-electron interaction operator, 
^ee'^ = Si<j ^ee''(^ij) is & loug-range electron- 
electron interaction operator with u'gg'''(r) = erf (/ir)/r, 
E^'^^ln] = (1/2) // n(ri)n(r2)wfe'^(ri2)dridr2 is a short- 
range Hartree functional with wf^'^{r) = erfc(/ir)/r, 
El'^i'ln] = ($[n]|H^|g^'^|$[n]) - E^'f'ln] is a short-range 
exchange functional where W|g'^ — X]i<j ''^ee'''(''y ) ^^'^ 
^[n] is the KS determinant, is a short-range cor- 

relation functional defined so that Eq. ([T|) is exact, and 
n>ii is the density coming from The minimizing wave 
function in Eq. ([T|) will be denoted by In these 



'Electronic address: ' savin@lct.jussieu.fr| 



equations, /i is a parameter controlling the range of the 
decomposition of the Coulomb interaction. In practice, 
approximations must be used for the wave function '^^'^'f^ 
and the short-range functionals E^'^''^ [n] and E^'^''^ [n] . In 
particular, the local density approximations have been 
constructed for El^'t'ln] and i?f'''[n] For /x = 0, 

Eq. ^ reduces to the KS scheme. In fact, in this case, 
the long-range interaction vanishes, M^ie'^^° = 0, 
the short-range functionals E^''^~ [n], El^'f^^^[n] and 
^sr,p=o jj^j j-g(j^gg ^}jg Hartree, exchange and corre- 
lation functionals of the KS theory. 

Former experience with Eq. ([T]) has shown that in gen- 
eral the quality of the wave-function ^f'''^ obtained with 
a given approximate functional is much better than that 
of the functional itself. To extract the maximum infor- 
mation from v]/i''.A'^ -y^re propose in this work to compute 
the ground-state energy as 

E = (M-'-^'^lf + Vrre + Weel^f''^^ + 4"mdK"-]> (2) 

where Wee = J2i<j ^I'^ij the full Coulomb interaction 
operator and Efe m^x^ is a new short-range correlation 
functional defined such as Eq. ^ used with the exact 
wave function vJ/I'^a* is exact. 

We note that Eq. ^ can formally be made self- 
consistent by generalizing the "optimized effective po- 
tential" (OEP) approach (see, e.g., Refs. El HI) to the 
multi-determinantal extension of the KS scheme 

E = inf {(M-'-^-^'Hlf + + W^ee|^"-'^H) + E'^^.^^A^v]} , 

(3) 

where the infinimum is search over one-electron poten- 
tials v{r), and ^'^'"•''[u] and Uy are, respectively, the 
ground-state multi-determinantal wave function and den- 
sity of f + Wll''' + J2iV{ri). If it exists, the mini- 
mizing potential is v{r) = w„e(r) -t- 5E^''^[n]/5n{r) + 
SE^J'f'lnl/Snir) + 5El''^'[n]/5n{r), establishing the link 
with Eq. ([1]). Eq. ([3]) reduces to the OEP formulation of 
the KS theory when /i = 0. In this case, in fact, ^^''^"=0 
is the KS determinant and the functional E'^^^^'^in] re- 
duces to the correlation functional of the KS theory. In 
practice, because the potential corresponding to the func- 
tional derivative of E'^'^^^n] is small, we expect that the 
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corrections on ^^'''f^ brought by self-consistency are negli- 
gible with respect to the errors due to the approximations 

This work is devoted to the study of the functional 
-^c'nidW' turn our attention entirely to the correla- 
tion energy, without combining it with an approximate 
functional for exchange. It is thus a pleasure to dedicate 
this paper to Professor Hermann StoU who has been a pi- 
oneer in the study and app lication of correlation energy 
density functional [13, ITsf . 

The paper is organized as follows. In Sec HH we dis- 
cuss the short-range correlation functional ^^'mdl"'] ^^'^ 



its relation to the functional E^'^'^[n]. In Sec. IIIII we 
construct a local density approximation for E^^^^ [n] . In 
Sec. IIVI we assess the accuracy of this approximation for 
the He atom. Sec. |V] contains our conclusions. Atomic 
units (a.u.) are used throughout this work. 



II. THE SHORT-RANGE CORRELATION 



FUNCTIONAL E 



c,md L 



The short-range correlation functional ^-^'^mdM 
Eq. ^ is defined with respect to the multi-determinantal 
wave fimction '5''^''', in contrast to the short-range corre- 
lation functional £'^'''^[n] in Eq. ([T]) defined with respect 
to the one-determinant wave function It is easy to 
see that these two functionals are related to each other 
through 



(4) 



where 



- ((^'■^■''M|T4^-'^|*"-'^[n]) - {<f[n]\W:i''^mn])) . (5) 

The quantity A''"''"''^ vanishes for n = and fj, — > oo. 

It is interesting to study the behavior of E^'^^^ [n] in the 
limit of a very short-range interaction, i.e. when /i — > oo. 
In this limit, the short-range interaction behaves as [13] 



(6) 



lea ding to the following asymptotic expansion of E^'^'^' [8i, 
MM 



El 



(7) 



where n2,c(r,r) is the correlation on-top pair density. 
The asymptotic expansion of A'''"''''"'^ as /j, — > c» is ob- 
tained similarly from its definition, Eq. leading to 



^Ir— sr,/i — >oo 



-^|n,,(r,r)dr + 0(-i 



(8) 



The first terms in Eqs. ^ and ^ cancel, and therefore 
-^Tmd decays at least as 1 //i"^ when /i ^ oo . 



III. LOCAL DENSITY APPROXIMATION 

A local density approximation (LDA) can be con- 
structed for El'^''^ 



c,md L 



-^c,md,LDAl"'J 



n[r)e 



c.mdjUnif 



.f{n{r))dr, 



(9) 



where the corresponding correlation energy per particle 
in the uniform electron gas e^'^j uniii''^) given by 



^c,md,unif("') 



= Cunif(") 



(10) 



In Eq. (fTO|) . e^'^unif(rt) is the correlation energy per par- 
ticle defining the LDA approximation for i?^'''''[ri] (see 
Refs. Isl [isl). and A|[j^f ''''' (n) is given by 



unit 



in) 



^c.unif 



ir,n)wll'''{r)ATTr^dr, (11) 



where 5c i(5iif(?', J^) is the correlation pair-distribution 
function of a uniform electron gas with long-range inter- 
action ?i;gg''(r) and density n. The correlation hole of this 
"long-range" electron gas is then given by n g]^'^nif{r, n) . 

Since an estimate of the energy e^'^f ("■) from coupled- 
cluster calculations is available [1, [l3| , we only need to 
compute the term ^^^if'^'^in) to build the LDA func- 
tional of Eq. ([9]). In order to estimate A[[j^f'^'^(n) we pro- 
ceed as follows. We first notice that wll''^{r) = eric(iJLr)/r 
in the integrand of Eq. ifTTj) only samples the part of 
.g^'^^(^if (r, n) corresponding to r < 1/^. For the standard 
uniform electron gas (with full interaction 1/r) the "ex- 
tended Overhauser model" [1^ proved to be able to yield 
accurate results for gc,unif('', n) in the short-range region 
defined by r < rs, where Vs — (47r n/3)^^/'^. We can thus 
use this simple model to calculate g^'^^^ifir, n) and to pro- 
duce an estimate for A[[j^f''^(ri) that should be reliable 
for /X- values for which /ir^ > 1. 

The scattering equations of the "extended Overhauser 
model" are widely explained in Refs. [lO, HH. Here we 
solved the same equations with the electron-electron in- 
teraction erf (/ir)/r screened by a sphere of radius of 
uniform positive charge density n and attracting the elec- 
trons with the same modified interaction. 



erf (/ir) 



^erf(Mlr--rl)^^,^ 



r' \<.rs 



(12) 

This potential is reported in the Appendix of Ref. ,22l . 
where it has been used for two-electron atoms with 
very accurate results for the corresponding short-range 
correlation energy. Vcs [r, rs , /i) is a screened potential 
that tends to the "Overhauser potential" [13, [1^ when 
/i — > oo, and which goes to zero when — > 0. As in the 
original work of Overhauser [2^ , the idea behind Eq. (fT^ 
is that the radius of the screening "hole" is exactly equal 
to r,. 
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1: A sample of the pair-correlation functions 3^'iuiif(^' 



of a uniform electron gas of density 3/(47rrg) (here = 2) 
with long-range-only interaction erf (/ir)/r obtained from the 
"extended Overhauser model" [lO] [see Eq. (|12p ]. The corre- 
lation hole is given by n g^'^^^f(r,n). 



A sample of the pair-correlation functions g^c unif{r, n) 
that we have obtained is reported in Fig. [1] As long 
as /i is not large we clearly see the absence of the cusp 
[d'c unif ~ ^) ^'^^ zero for a system with interaction 
1 /r at small r , but it is zero for the erf (/xr) /r interaction] . 
As expected, as fi increases the hole deepens, and for very 
large /i we see that the cusp starts to appear. 

Some of the values of A^j^f'^'''(n) for fiVg > 1 are re- 
ported in Table HI An estimate of AJ^^^j'^''^ (n) in the re- 
gion not accessible with the extended Overhauser model, 
f^''"s ^ 1: has been obtained by a simple interpolation be- 
tween our data and zero, since, as explained in Sec. [TTl 
^ir-sr,p vanishes when fi ^ 0. In the opposite limit. 



system of uniform density reads 



00, A"j^|f'^'''(n) behaves as in Eq. which for a 



A: 



Ir— sr,/j-- 
unif 



8 fi'^ 



O 



(13) 



where (7c,unif(0, n) is the on-top value (r = 0) of the pair- 
correlation function of the Coulombic uniform electron 
gas of density n We found that the A[[j^j'^'''(7i) com- 
puted with the extended Overhauser model accurately 
recover this limiting behavior. 

For future applications, a more accurate LDA func- 
tional for A^''"'^''''^ (especially for iiVg ^ 1) will be avail- 
able from quantum Monte Carlo calculations [23| . 



IV. RESULTS FOR THE HE ATOM 

For the He atom, the short-range correlation energies 
E"^''^ and E^'^j^^ have been calculated with a precision 
of the order of 1 mH as follows. An accurate density is 
calculated at the full configuration interaction level with 



a large Gaussian basis set and the optimization of the 
potential in the Legendre transform formulation (25l . [26j 
of density functionals enables to compute accurately the 
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FIG. 2: Short-range correlation energies i?^"^'*^ (thin curves) 
and (thick curves) with respect to /j, for the He atom. 

Accurate calculations (solid curves) are compared to the LDA 
approximation (dashed curves). 



correlation energy E^'^''^ associated to that density (see 
Refs. m, [3 [13 for details). The corresponding accurate 
multi-determinantal wave function \I>''''^ and the KS wave 
function $ are also obtained in this procedure, which give 
access to an accurate evaluation of A'"'"*'''''' and conse- 
quently of E^^^:^. 

The accurate and LDA short-range correlation energies 
E^'^'^^ and i^^'^i^d ^'^^ compared in Fig. [2l For all values 
of /i, we have |i?^'mdl < l^f^^L meaning that in Eq. © 
(^k,M|iy|r,;.|^ir,A.^ < {^iWll-t^l^) which secms natural 
for a repulsive interaction. Tab. |TT] compare the LDA 
errors on E^^'^ and i?^'md- ^^ne sees that the LDA errors 
for this two short-range correlation energies are of the 
same order of magnitude for all values of /x. 



V. CONCLUSIONS 

In this work, we have reexamined the multi- 
determinantal extension of the KS scheme based on a 
long-range/short-range decomposition of the Coulomb 
electron-electron interaction. Contrary to previous works 
where the short-range correlation functional was defined 
with respect to the KS determinant, we have introduced 
a new short-range correlation functional defined with re- 
spect to the multi-determinantal wave function. We have 
constructed the local density approximation for this new 
functional. The example of the He atom suggests that 
the local density approximation is essentially as accurate 
as for the short-range correlation functional defined with 
respect to the KS determinant. We believe that this work 
paves the way to a multi-determinantal extension of the 
KS scheme using a correlation-only density functional. 
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M ^unif 


M '^s ^unif 


M ^unif 


A Ir — sr.Lt 


II r Alr-sr,M 
M ^unif 


2 0.5 0.0173 
2 1 0.00979 
2 2 0.00310 
2 3 0.00125 
2 4 0.000608 
2 5 0.000335 
2 6 0.000202 


3 0.5 0.0136 
3 1 0.00585 
3 2 0.00153 
3 3 0.000586 
3 4 0.000278 
3 5 0.000151 
3 6 9.08 ■ 10-^ 


5 0.2 0.0175 
5 0.5 0.00813 
5 1 0.00265 
5 2 0.000598 
5 3 0.000219 
5 4 0.000102 
5 5 5.52 • 10^^ 


10 0.2 0.0109 
10 0.5 0.00299 
10 1 0.000776 
10 2 0.000158 
10 3 5.63 • 10"^ 
10 4 2.59 ■ 10"^ 
10 5 1.39 ■ IQ-^ 


15 0.2 0.00683 
15 0.5 0.00151 
15 1 0.000363 
15 2 7.16 ■ lO"'^ 
15 3 2.52 • 10"^ 
15 4 1.16 ■ 10"^ 
15 5 6.19 ■ 10"^ 



TABLE I: A sample of the values of A[[^;j'^''^ [see Eqs. (|10|) - (|1H )] computed from the extended Overhauser model [20| . 
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TABLE IL LDA errors AE'f ''^ = S^'^a ^ Ef" and AEl';^^^ = K'Xd.i^uA ' K]u^d with respect to ^ for the He atom. 
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